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Abstract 

We consider transformations of the 2x2 propagator matrix in 
real-time finite-temperature field theory, resulting in transformed n- 
point functions. As special cases of such a transformation we examine 
the Keldysh basis, the retarded/advanced RA basis, and a Feynman- 
like FF basis, which differ in this context as to how "economically" 
certain constraints on the original propagator matrix elements are 







implemented. We also obtain the relation between some of these real- 
time functions and certain analytic continuations of the imaginary- 
time functions. Finally, we compare some aspects of these bases which 
arise in practical calculations. 
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1 Introduction 



The fields = x) describing particles in a heat bath at temperature 
f3~^ are either periodic (bosons) or anti-periodic (fermions) in the complex 
time plane, with period —if] [1-5]. Integrals over time then follow a path C 
in the complex time plane connecting a point — T to a point — T — i(3. In 
the imaginary-time formalism the straight path connecting these points is 
chosen, resulting in time arguments t = —it that are purely imaginary - an 
analytic continuation is then involved for quantities defined for real times. 

In real-time formalisms a path which includes the real time axis is chosen, 
and thermal Green functions are then time-ordered with respect to this path: 

Gc{xi, . . . , x„) = (Tc0(xi) ■ ■ ■ (t){xn)) . (1) 

A standard choice of this path is the concatenation of the following four 
pieces: Cj : —T^T; Cm: T^T — icr; Cu '■ T — ia ^ —T — ia; Cjv '■ 
—T — ia — T — if3, where a is an arbitrary parameter < a < /? 0, ^]. 
As T ^ oo the contributions from Cui and C/y can be neglected in the 
majority of cases. It is then convenient to introduce the thermal doublet 

where a,b = +,—■ The fields (p+{x) = and 4>-{x) = (pit — ia,x) are 
defined on the forward and backward segments of the time contour, Cj and 
Cu respectively. This leads to a matrix form of the thermal Green function: 

Ca^...a,Xxi, . . . ,X„) = (Tc0ai(Xl) ■ ■ ■ </>a„(a;„)) . (3) 

The dependence on a will not be indicated explicitly, except in cases where 
a particular value applies. Two popular choices for this contour parameter 
that have been studied are a = 0, which is the closed time-path formalism of 
Schwinger, Keldysh and others [6-9], and a = \, which is equivalent to the 
"unitary" formulation of thermo field dynamics ||^. 

The propagator for free particles is a 2 x 2 matrix and in momentum 
space it has the form 

= f ^++(^^) (4) 
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where 

= D*__ip) = 

= e{po)AF{p) - e{~po)A*M + veiPo)n{x)[AF{p) + A*^, 

= r)e{po)n{x)e^P^^[Ap{p) + A*p{p)]. (5) 

Here Af{p) is the Feynman propagator — + ie] for scalar bosons), 

and ^ 

is the Bose-Einstein (Fermi-Dirac) distribution function with the properties 

e^n{x) = —rin{—x), 

n{x) + n{—x) = —1], (7) 

where x = P{po — /i is the chemical potential, and t] = +1 (—1) for bosons 
(fermions) . 

In many applications Green functions other than these time-ordered prod- 
ucts arise. A common example is the linear response of the field to a weak 
external current ^, which involves the retarded propagator A/j. This 
product also arises in a comparison of certain analytic continuations of the 
ra-point functions in the imaginary-time formalism to functions in real-time 
formalisms [|, [Tl|, [l^, |13[. As the retarded propagator can be written as 

An{p) = D++ip) - e-^^'D+4p) = e(po)A^(p) - 9{-po)AUp), (8) 

which for bosons is i/lp"^ — + ispo]-, it is then of interest to consider trans- 
formations of Eq. @) which could generate such linear combinations of the 

Dab. _ 

In this paper we will examine transformations to a Dxy basis of the form 

U{p)D{p)U^{-p) = D{p), (9) 
where "T" denotes the transpose and 
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Without further input, however, there would hkely be no advantage to study 
such a transformation in generaL One observation that could be a guide 
for choosing "useful" transformations is that since the original propagators 
Dabij)) satisfy the relations 

+ Z}__(p) - e-'^P«D+_(p) - e'^P«Z}_+(p) = 0, (11) 
D++{p) + D__{p) - r7e'^P«-"D_+(p) - T^e^'-'^^W+^p) = 0, (12) 

the transformed Dxy{p) of Eq- (H) will also not be independent in general. 
It could then be of some advantage to choose the transformation so that 
either one or both of the relations of Eqs. (0, |12D translate into the trivial 
vanishing of either one or two of the Dxy{p)- 

It may be noted that the transformation of Eq. (^ has a close correspon- 
dence to the Bogoliubov transformation studied in thermo field dynamics 



14]. In fact, the transformation matrices U{p) can be mapped onto a par- 
ticular set of Bogoliubov transformations seeking a diagonal form for the 
propagator by a transformation of the field operators p| . 



2 Transformed n— point functions 

The first relation of Eq. ( [Tl| ) among the Dab{p) elements is an immediate 
consequence of the definition of the path-ordered product, whereas for the 
second one of Eq. (|12|) one also needs the Kubo-Martin-Schwinger (KMS) 
condition for equilibrium two-point correlation functions. These relations are 
thus also satisfied by the full propagator Gabip)- Hence, the full propagator 
has the same matrix components of Eq. (^ as the bare one, with now the 
Feynman propagators in Eq. (^ given by the spectral representation 

JO Pq — + ie 
This implies that for the full propagator we may write the spectral form 

/•oo 

Gabip) = dT^p{T,p)Dab{T,p) (14) 

Jo 

and the analogous one for the transformed propagator 

^ POO , 

Gxy{p)=\ dT^piT,p)DxYiT,p), (15) 

Jo 
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owing to the identical matrix structures of interacting and non-interacting 
propagators. 

Consider now Dyson's equation for the full propagator defined on the 
contour C/ U C// g |]: 

/ dyTc{x,y)Gc{y,z) = 6c{x - z), (16) 

or, in a matrix representation, 

/ + 00 
dyT{x,y)a3G{y,z) = a36{x - z), (17) 
-oo 

where cxs = diag(l, — 1) is the 2x2 Pauli matrix. Writing Tc{x,y) = 
r[?^(x, y) — Tic{x, y), in momentum space we then have 

Gip) = Dip) + Dip)J:ip)G{p), (18) 

where here, and also in the following for n-point functions, the effects of 
have been absorbed into the vertex function by associating a minus sign with 
each "— " index. From Eq. (P) and the preceding discussion we thus see 

U{p)G{p)U^{-p) = G{p), (19) 

where G{p) satisfies 

Gip) = Dip) + Dip)±ip)Gip) (20) 
and the transformed (amputated) self-energy S(p) is given by 

v{p)j:ip)v^i-p) = s(p), (21) 

where we have defined 

^""^^^^ - \ v,^{p) v,.{p) )-[^ ^ m ■ ^^^J 

Generalizing this we can, as with Aurenche and Becherrawy |]15|, |l^], 
introduce transformed (amputated) n-point functions V x^-x„{pi, ■ ■ ■ ,Pn) in 
terms of the +/— functions Taj^...a„{pi, . . . ,Pn)- If all momenta are incoming, 
such functions are 

fxi xAPU ...,Pn) = Vx.aiiPl) " " " Vx^aniPn)'^ ai-aAPl^ ■ ■ ■ ^ Pn) ■ (23) 
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If some of the momenta were outgoing, say Pk+i, ■ ■ ■ ,Pn, the matrix V^l—p) 
rather than V{p) would be associated with these momenta in Eq. (|23|) . Be- 
cause of this, the corresponding transformed n-point function is then simply 
related to that with all momenta incoming by 

fxi.-Xfe;Xfe+i-..X„(Pl, • • • ,Pfc;Pfc+l, . . . ,Pn) = Txi-X„iPl, ■ ■ ■ ,Pk, -Pk+1, • • • , -Pn 

(24) 

Thus, it suffices to study these functions with all momenta incoming. We 
will assume this in the following, so that J2Pio = = and we also 

use fermion number conservation: Ylrji = 1. Also in what follows the delta 
functions associated with these conservation laws have been factored out 
from the Green functions. 

Feynman rules for the transformed functions of Eq. (p3D then follow from 
the rules for the n-point functions ra^...a^, which involve internal propagators 
Dab connected by appropriate bare vertices gab-'- 

g+...+ {pi, ...,pn) = -g-...-{pi, ...,pn)= g{pu ■ ■ .,Pn), (25) 

with all others being zero [0. Instead of these vertices one may also consider 
amputated subdiagrams, for which the following discussion will also apply. 
The internal propagators in a graph could be written in terms of transformed 
propagators Dxy using Eq. @ and, as also done by Aurenche and Becher- 
rawy [|T3|, the matrices of each such transformation could be associated with 
the adjacent vertices. This will ultimately replace each vertex gab -- in a 
graph by the appropriate transformed vertex gxY---, where the transforma- 
tion is given in Eq. (pSj). The end result is that the Feynman rules for the 
transformed n-point functions can be written solely in terms of transformed 
vertices connected by transformed internal propagators. 

A simple illustration of this transformation of Feynman rules is a one-loop 
self-energy graph involving cubic interactions such as 

f d^k 

-iEaa'ip) = (-0 y gabcip, -k,q)Dbb'{k)ga'b'c'{-p,k, -q)Dc'c{q), 

(26) 

where q = k — p and only thermal indices are retained. The transformed 
self-energy S(p) of Eq. ( ^1]) corresponding to this graph is, using Eq. (^, 

-iExX'ip) = Vxa{p) [-iT.aa'{p)] Vx'a'{-p) = 

= {-if J gxYz{p, -k, q)DYY'{k)gx'Y'z'{-p, k, -q)Dz'z{.<W) 
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where, as in Eq. (|23|), 



gxYzip, -k, q) = Vxa{p)VYb{-k)Vzc{q)gabc{p, -k, q), 

9X'Y'Z'{-P, k, -q) = Vx>a'{-p)VY>b'{k)Vz'c'{-q)ga'b'c'{-p, k, -g).(28) 

It will be found convenient in the construction of the following transfor- 
mations to aim for the +/— n-point functions rai...a„ to be multiplied by 
one of two particular factors for each index with the value "— " : 



Ta,...aAPl,---,Pn) H 

a-i=— 

Ta,...a,Xpl,...,Pn) Y[ T^.e'^^'O "^^ • 



(29) 



The exponential factors exactly eliminate the cr-dependence of raj...a„, leav- 
ing these combinations a-independent. The first function here is simply the 
(7 = n-point function expressed in terms of Taj^...a„ with an arbitrary a. 
The second function can be related to the first by the KMS relation. These 
functions satisfy 



E 

ai=± 

E 

ai=± 



Ta,...aAPU---,Pn) H = 0' 

ai=- 

Ta,...a,APU---,Pn) U //.c'^^'O-- 



(30) 



which can be proven, for example, using the largest and smallest time equa- 



tions of 't Hooft and Veltman |T3, |18|. These relations are the amputated 
versions of the n-point generalizations of the 2-point relations of Eqs. (|lT], 
|l^). Note that, unlike Eqs. ( PI |12D , in Eq. ( ^OD no minus signs appear due to 
the their absorption into the amputated vertex functions, as discussed after 
Eq. (0). 

We shall also have occasion to consider the particular combinations 



r(^)(pi;P2,...,Pn)= E 



ai=± 



r(^)(pi;P2,...,Pn)= E 

ai=± 



r+,,...,„(pi,...,p„) n e''^'' 

ai=— 



f3ll 
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which are real-time representations of the one-particle-irreducible contribu- 
tions to the generalized n-point retarded and advanced products |16|, |l^]. 
They are related by complex conjugation: T^"^\pi;p2, . . . ,Pn) = T^^-* *{Pi',P2, , • • • ,Pn}, 
which can be shown using Eq. (^) together with the conjugation relation 



r:,....„(Pi, . . . ,Pn) = -r.,....„(pi, . . . n ^.e^-^'^^'o, (32) 

ai=- 

where a = — , + when a = +, As will be seen later, the identification of 
r*^-^^ and r*^"^) with the retarded and advanced amputated functions, respec- 
tively, is based on the convention in the imaginary-time formalism for the 
retarded and advanced analytic continuations of the (non-amputated) Green 



functions [lIOl 0, |l2l IT 



3 The Keldysh basis 

The Keldysh transformation of the time-path method is an example where 
one of the relations of Eqs. (0, |l^) is used to eliminate one of the trans- 
formed functions Dxy{p) [6^9]. Such a transformation is particularly useful 
in non-equilibrium situations, where a relation such as Eq. (p!2[) , which fol- 
lows from the KMS condition, does not hold. We will find that this class of 
transformations automatically involves the retarded and advanced propaga- 
tors A/j(p) = —A\{p) of Eq. (^, but there is some freedom in the form of 
the third non-zero element. As is conventional, we choose this element to be 

proportional to the combination D^^{p) + D (p). 

We consider the case where the transformed function D22{p) of Eq. (|^) 
vanishes identically by Eq. (|ll]). With fcj = 1, 2, the transformed propaga- 
tors, 

with the normalization on the right-hand-side, follow by choosing in Eq. (|^) 

_ 1 / b{p)b{-p) e'^P<^b{p)b{-p) \ 
'"^^^ ~ V2b{-p) [ 1 -e'^P" ) ' ^^^^ 

where b{p) is an arbitrary function and the symmetric product As{p) is 

As{p) = b{p)b{-p) [D++{p) + I^— (p)] = b{p)b{-p) coth''(a;/2) [Ar{p) - Aa{p)] 

(35) 
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We then define "1/2" n-point functions ICki-.-kSvii ■ ■ ■ iVn) in terms of the 
+/— n-point functions Va^...a^{j)i, . . . ,pn) by Eq. (p3l). The corresponding 
matrix V{p) of Eq. in this case is, from Eq. (P^, 

/ \ ^ 1 / 1 e-^° \ 

1^ V2+(p) V^2-(P) y V2h{-p) \ b{p)b{-p) -e^v^h{p)h{-p) ) ■ ^^"^ 

In this notation the amputated 2-point functions are 

/Cll(pi,P2) = 0, 

^12(^1,^2) = r(^)(pi;p2), 

/C2i(pi,P2) = r(^npi;P2), 

/C22(pi,P2) = 6(pi)6(P2)coth''(xp,/2) [r(«)(pi;p2) - r(^)(pi;p2)] ,(37) 

and in general one can show, using the largest /smallest time equation rela- 
tions of Eq. (|30|) , that 

/Cil...l(pi, . . .,Pn) = 0, 

/C211...1(P1, . . . ,Pn) = 2^-/2 ^^^4^ rr(^nPi;P2, . . . ,Pn), (38) 

6(-p2)---K-Pn) 

where T^^\pi;p2, . . . ,Pn) is defined in Eq. (pT]). 

Feynman rules for the 1/2 functions consist of internal propagators of 
Eq. (^) and bare vertices obtained from Eqs. ( p3| , |36D using the bare +/— 
vertices of Eq. (ESl) : 



where "#2" denotes the number of "2" indices present. Thus, bare vertices 
with an even number of "2" indices vanish. The non-vanishing ones satisfy 
the following relation when the signs of all energies and momenta are reversed: 

9kr-kr,{Pl,---,Pn) g{pi,...,Pn) ^^[/[^^ ^Pi) 

In these considerations there is no obvious advantage to choosing the free 
parameter b{p) to be momentum dependent, and so choices such as b{p) = 1 
or I/V2 are usually most convenient. 
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As in this type of transformation only the one constraint of Eq. (^) 
satisfied by the Dab was used to ehminate one of the transformed Green 
functions, the other constraint of Eq. ([T2|) will lead to non-trivial relations 
amongst the transformed functions Kk^...^^. This is illustrated by the last 
relation of Eq. (|37|) relating /C22 to /C12 and /C2I) which is referred to as the 
fluctuation-dissipation theorem 0. 

One could also consider the case where the transformed function -D22(p) 
of Eq. (P) vanishes by Eq. (0), rather than by Eq. (|ll|) as chosen here. The 
analysis of such a transformation is straightforward; we only note that in this 
case the transformed n~point function /C2ii - i will be proportional to V'^'^\ 
rather than to V^^^ as happened in Eq. (^). 

These Keldysh-like transformations are convenient in a general analysis 
of response theory |^ . In this approach one adds to the Hamiltonian a term 
involving a (weak) external current J{t) and then calculates the response of 
various Green functions. Translating from the notation of Ref. into that 
used here, one has 

5Gi{t)r^ J dh Gu{t,h)J{ti) + ^J dtidt2 Gu2{t,h,t2)J{h)J{t2) + ..., 
SGn{t,t')^ J dtiGn2{t,t',ti)J{h) + ..., (41) 

and so on. The Green functions G'i22- -2(^i, ■ ■ ■ ,tn) with one "1" index are the 
n-point retarded products, and arise in the response of the expectation value 
Gi (t) of the field operator, while those response functions with more than one 
"1" index arise in the response of higher-order correlation functions G11...1. 
The one-particle irreducible functions )Cki---kn considered in this section enter 
into the construction of these response functions, and as such this setting 
provides an intuitive interpretation of them. 

4 The RA basis 

In equilibrium situations one can consider transformations in which both of 
the constraints of Eqs. ([TT|, p^) translate into the trivial vanishing of two of 
the transformed functions Dxy{p) of Eq. (P). This might be viewed as an 
attempt to "economize" the information contained in the original propaga- 
tors Dab{p)- We will consider two classes of such transformations. We first 
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look at the one that results by demanding Du(p) vanishes by Eq. (|Il|) and 
D22{p) vanishes by Eq. ([T^). We find in this case that the transformation 
automatically leads to an off-diagonalization in terms of the retarded and 



advanced propagators |16 



where ai = R, A and, for the normalization on the right-hand-side to occur, 

.X _ 1 / a{p)a{-p) -e'^Poa(p)a(-p) \ . . 

with a{p) an arbitrary function. 

By Eq. ( ^31) we then define ''RA" n-point functions 7^q,^...q,„(pi, . . . ,Pn)] 
the matrix V{p) of Eq. ( p2D is, from Eq. (^), 



^R+ip) ^R-ip) \ _ 1 / —r]n{—x) 7]e"'P°n{x) 
VA+ip) VA-ip) j " <^ V e-foa(p)a(-p) 



(44) 



In this notation the non-zero 2-point functions are the self-energies 'JZra{pi,P2) 
TZ\^{pi,P2) = T^^\pi;p2). Using the largest /smallest time equation relations 
of Eq. (0), one finds in general [|l^, 

T^RR-RiPl, • • • ,Pn) = = 71aA-a{Pi, ■ ■ ■ , Pn) , 

I^-RAA -AKPi-, . . . ,P„) = 7 \P\\P2-, ■ ■ ■ ,Pn), 

a{-pi) 

nARR...R{Pl,---,Pn)-{-l) n{x,M-P2)---a{-Pn) ' ' " " " ' ' 

(45) 

where T'^^\pi;p2, . . . ,Pn) and T'^^\pi;p2, . . . ,Pn) appear in Eq. (|3TD. 

The RA functions satisfy a complex conjugation relation; relating Vaa{p) 
to Vaa{p) and using Eq. (|32D, one finds 



Ki...a„(Pl,P2, . . . ,P„) = -7^al...a„(Pl,P2, • • • ,Pn) X 

-pj- n(x^) -j-j- a(p»)a(-p^) 
}=R a{pi)a{-pi) ^^i^ -n(-Xi) 
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where a = A, R for a = R, A. This agrees with the results of Refs. |]T3|, |T6| 
for the choice a{p) = —rjn{x). We see by Eq. (^6[) that it is not possible 
to choose a{p) so as to cancel the factors of the distribution function n{x) 
arising in this complex conjugation relation. 

Feynman rules for the RA functions consist of internal retarded or ad- 
vanced propagators and bare vertices obtained from Eqs. (p3| , |4^ ) using the 
bare +/— vertices of Eq. (^Sf ): 



9ai-a„{Pl, ■■■,Pn) 



9{P 



1) 



,Pn) 



n 



R a[-Pt) a,=A 



n 



n - n 



ai=R 



ai=R 



(47) 



Thus, as for the full ra-point functions in Eq. (^), bare vertices with all "i?" 
or all "A" indices vanish: gRR...R = = gAA---A- The non-vanishing ones 
satisfy 



-Pu 



-Pn 



9ai-aAPU 



,Pn 



-Pnj 



g{pi,...,Pn) 



TT O'i.-Pi) 
iL a{pi) 



(48) 



where ^R denotes the number of "i?" indices present. 

Two convenient choices of the free parameter a{p) can be made based 
on these considerations. One is a{p) = 1, which normalizes the bare ver- 
tex gRAA-- A to the bare +/— vertex g of Eq. (|25|). In this case the factor 
Ila{—p)/a{p) in Eq. (|48|) is unity, and the full n-point function 71raa--a of 
Eq. (1^) is normalized to T^^\ The other choice, as made in Refs. [l^, [l^] , 
is a{p) = —rin{x), which normalizes (Jarr - r. to the bare +/— vertex g. For 
this choice the factor I\a{—p)/a{p) in Eq. (^Sj) is (—1)"", and the full n-point 
function TZarr-r of Eq. (|^) is normalized to F*^"^). 



5 The FF basis 

As a second example of using Eqs. (0, |T2|) to eliminate two of the transformed 
functions Dxrip) in Eq- (i), "we consider the case where the off-diagonal 
components D^lp) and D2i{p) vanish. To avoid a singular transformation 
we find we have to impose, for example, that D2i{p) vanishes by Eq. (|ll]) 
for po > and by Eq. (|I2[) for po < 0, and vice-versa for Di2{p). Such 
a transformation automatically leads to a diagonalization in terms of the 
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Feynman and anti-Feynman propagators: 



DhfM 



Dp rip) Dp pip) 
Dpp{p) Dpp{p) 



Ap{p) 
AUp) 



(49) 



where fi = F,F and, for the normahzation on the right-hand-side to occur, 

Ufaip) = eipo) 



c{p) — 776°"^° ^c{p) 
-rjdijp) rje'''P°d{p) 



n[ — X e 



-O-PO 



d{-p) -d{-p) 



c{-p)d{-p) V -e"-'^P«c(-p) r]c{-p) j ' 

, , pn) by Eq. 



(50) 



where c{p) and d{p) are arbitrary functions. 

We next introduce "FF" n-point functions ^/i - /„(pi 
the matrix V{p) of Eq. (p2D is, from Eq. (0), 



VpM Vp^{p) 

VpM yp-ip) 



0{Po 



c{p) rie"^° ''^c{p) 
r]d{p) r]e^^°d{p) 



+e{-po) 



n(—x)e 



apo —X 



e~''P'>d{-p) d{-p) 

c{-p)d{-p) \ e-"^P<^ci-p) r]ci-p) 



(51) 



The non-vanishing 2-point functions are then J-'pp{pi,P2) = —^p'p{Pi,P2) = 
9{pif^)T^^\pi; P2) + 9{—pig)T^^\pi; P2) . It is worth noting that J-'pp, not the 
original amputated 2-point function r_(_+, determines the pole of the time- 
ordered Green function G++(p) of Eq. (|^), and that at finite temperature 
these two functions do not coincide, even at the one-loop level [|, 0]. A 
similar lack of equivalence can also be shown for the higher n-point functions 

J'pp...p and r_|_^ |_. 

Feynman rules for the FF functions consist of the Feynman and anti- 
Feynman propagators and bare FF vertices obtained from Eqs. (|2^, ^) 
using the bare +/— vertices of Eq. (^S]): 



dh-u{Pi,---,Pn)/g{pi,- 
9{Pio)c{pi) + B{-pi^ 



n 
n 



0{Pio)c{Pi) + Vi9{-Pi, 



■,Pn) = 

c{-Pi) 
n{-Xi) 



n 



C{-Pt 



n 

J Si=F L 



m^iPiMiPi) + ^{-Pi. 

0{Pio)e^'d{pi) + 9{-pi^ 



'1-0 ) J/ 

d{-Pi, 
d{-pi) _ 



(52) 
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which satisfy 



9h-fA-Pi 



•I • • • -I Pn 



n{xi) 



n 

f^=F,F L 



OiPn 



gj-pu - ■ ■ ,-Pn) ^ 
g{pi,...,Pn) 

c{-pi)d{-pi) 



c{pi)d{pi) 



+ 0{-p,,)- 



n{-Xi) 



(53) 



Two considerations illustrate the relative effects of different choices of the 
arbitrary coefficients c{p) and d{p) in this basis. One is the normalization of 
the FF n-point functions, for which we have, for example. 



0{h)e{-k) ■ ■ ■ 9i-ro)J'FF..Ak, r) = 

c{k)n{—xi) ■ ■ ■ n{—Xr) 
_ n{xk)c{-l) ■ ■ ■c(-r) 
0{ko)ei-lo) ■ ■ ■ ei-qo)6{ro)J^pF...Mk, l,...,q,r) 



9iko)ei-lo)---e{-ro) 



r(«)(A;;/,...,r), 



e(ko)9i-lo)---ei-qo)eiro) 



c{k)n{—xi) ■ ■ ■n{—Xq)d{r) 



_n{xk)c{-l) ■ ■ ■ c{~q)rire-''- 
'd{l)---d{ry 



-e{^ko)e{k)---e{ro) , 

-ko)e{lo) ■ ■ ■ e{ro)TFF-F{k, l,...,r) 
c(/) ■ ■ ■ c(r) 



r(^)(A;;/,...,r), 



e{-ko)eiio)---eir,) 



ci-k) 

9{-ko)9ilo)---9{qo)9{-ro) 



9{-k,)9{k) ■ ■ ■ eiqo)9i-ro)TFp...Mk, l....,q,r)-- 

c(Z) ■ ■ ■ c{q)n{—Xr) 
c{—k)d{—r) 
9{ko)9{-k) ■ ■ ■ 9{-r^)J^pp...p{k, /,..., r) = 



-^^(A;o)^^(-/o)---^(-ro) 



d{k)n{—xi) ■ • ■ n{—Xr) 
n{xk)d{—l) ■ ■ ■ d{—r) 



T^^\k;l,...,q,r), 



r(^)(fc; /,..., r). 



(54) 



where T^^\k; I, . . . ,r) and T^^\k; I, . . . ,r) are defined in Eq. (|3lD. The other 
consideration is complex conjugation: relating Vja to Vf^ from Eq. (^T]) and 
using Eq. (^^ one can show 

^/V-/n(Pl' ■■■,Pn) = -^h-uiPU ■ . . ,Pn) X 
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n 



d{pi) 



ci-Pi) 



n 



d{-pi 



(55) 



where f = F, F for f = F, F. 

We consider now in the preceding context two choices of the parameters 
c{p) and d{p). For the first one, we begin by noting that the "symmetrical" 
or "unitary" formulation of thermo field dynamics is recovered by the choices 



a 



c{p) = e^/yn(a;), 
d{p) =r7e"^^/2yn(^ 



(56) 



in the transformation matrix of Eq. (^Tj) ^ Of course, potential 

problems exist for massive bosons in this when taking the square root of the 
distribution function. However, the 6 functions of energy which accompany 
c{p) and d{p) in general guarantee that this square root remains real except 
in the case of finite chemical potential with < {po/ fi) < 1, where a separate 
analysis in this and the following is required 0, |T5| . Guided by the values of 
Eq. (|56D we then consider, for arbitrary a, the choices 



c{p) 



n{x), 



d{p) = 7]^ n{x). 
In this case the bare vertices of Eq. (|5^) are 

9h-fSPi^---^Pn) 



9{Pi, ■ ■ 
which satisfy 

and where 



,Pn) 



n c[p^] n v^s[p^] - n ^b^i n c[p^ 



9h-fn{-Pl, 



-Pn 



9[-Pi, 



-Pn) 



9h-fAPl, - ■ ■ ,Pn) g{pi,...,pn) 

C[p] = e{p^)e^'^^^) + e{-p^)e--'^^n{-x) 
S[p] = e{po)^n{x) + rie{-po)\/n{-x), 



(57) 



(58) 



(59) 



(60) 



15 



which obey C^[p] — riS'^[p\ = 1. This choice of parameters leads to the 
following prefactors in the normalization relations of Eq. (^41): 



c{k)n{-xi) ■ 


■ ■ n 


{-Xr) 


n{-xi) ■ 


■ ■ n{ 


Xr) 


n{xk)c{-l) 


■ ■ ■ c 


i-r) ~\ 


n{xk) 


c{k)n{-xi) 


■ ■ n 


{~Xg)d{r) 




n{- 


-Xi) ■ ■ ■ n{—Xq)n{xr) 


n{xk)c{-T) ■ 


■•c( 




' \ 




n{xk) 


d{l) ■ ■ ■ d{r) 




n{xi) ■ ■ ■ n{xr) 






d{-k) 




n{-Xk) 






c{l) ■ ■ ■ c(r) 




n{xi) ■ ■ - n 


[Xr] 







c{-k) ^ n{-Xk) 

c{l) ■ ■ ■ c{q)n{-Xr) ^ ^ ^_^,,/2 
c{—k)d{—r) ^ 

d{k)n{—xi) ■ ■ -n^—Xr) 



n{xi) ■ ■ ■ n{xq)n{—Xr 
^ n{-Xk) 



n{—xi) ■ ■ ■ n{—Xr) 
n{xk)d{-l) ■ ■ ■ d{-r) ^ n{xk) 

while the complex conjugation relation of Eq. (|55|) is 



^h-fS^^ • • • = -^h-u(k, . . . ,r) n ^i- 

f^=P' 

As a second choice of the parameters c{p) and d{p), we consider 

c{p) = e^n{x), 
d{p) = 7]. 

This choice leads to bare vertices of Eq. ( |52D of the form 



(61) 



(62) 



(63) 



/ N I[ + Vi^iPioMxi)] Y[ [9{pi^) + rii0{-pi^)n{-Xi)] - 



n [vSi.Pia)'^i.^i) + ^i.-Pio)] n + ViG{.-PiM{,-Xi)] 

h=F fi=F 



(64) 



which satisfy 



9h-fA-Pi, 



-Pn) gi-Pi,---,-Pn) 



9h-fSPl^---^Pn) 



gipu---,Pv 



(65) 



16 



In this case the prefactors in the normahzation conditions of Eq. (|5^ ) become 

c{k)n{—Xi)---n{—Xr) ^ 
n{xk)c{-l) ■ ■ ■ c{-r) 
c{k)n{—Xi)---n{—Xq)d{r) 



n{xk)c{-l) ■ ■ ■c{-q)r]re 
d{l)---d{r) _ 

d{-k) ' 
c(/)---c(r) n{xi) ■ ■ -ni^Xr) 



c{-k) n{-Xk) 
c{l) ■ ■ ■ c{q)n{—Xr 



ji{xi) ■ ■ ■ n{xq)n{—Xr 



c{—k)d{—r) n{—Xk) 
d{k)n{—xi) ■ ■ ■ n{—Xr) n{—xi) ■ ■ ■ n{—Xr) 
n{xk)d{-l) ■ ■ ■ d{-r) n{xk) 

while the complex conjugation relation of Eq. (|55|) is 

•^/i-/n(Pl' • • • ,Pn) = -^h-fSVl. • • • ,Pn) X 



(66) 



n 



0{pio)n{xi) + 



n{-Xi) 



n 

fi=F L 



niXi 



+ 9{~pi^^n{-Xi 



One can also consider parameters c{p) and d{p) that represent choices 
that interpolate between the two competing features of normalization and 
complex conjugation. For example, the choice c(p) = e^/^ and d{p) = 1 
leads to a simple complex conjugation relation of Eq. ( pSj) but also to a 
somewhat more involved set of normalization relations of Eq. (|54D ; this has 
the advantage that no square roots of distribution functions need be taken, 
but at the expense of having a more complicated form of Eq. (|53D relating 
the bare vertex . . . , -Pn) to • • • ,Pn)- 



6 Comparisons of the bases 

In this section we look at the transformations between the three bases con- 
sidered here. As well as comparing their respective structures, this will allow 
us to see how relations in one basis translate into relations in another ba- 
sis. We first consider the transformation between the RA functions and the 
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Keldysh-like 1/2 functions. Writing this relation as 



n 



XPi, . . .,Pn) = Ta.kAPl) ■ ■■Ta,,k,XPn)}Ck,...k„{Pl^ ■ ■ ■ , Pn) (68) 



we find, for the choice b(p) = 1/ V2 and a{p) = 1 in Eqs. the matrix 

(69) 



Tri{p) Tr2{p) 
Tai{p) Ta2{,p) 



|coth"(s/2) 1 
1 



With this transformation, it is seen that the identity TZaa - a = in the RA 
basis translates in the Keldysh basis to /C11...1 = 0. The other identity in the 
RA basis, TZrr...r = 0, gives a relation amongst various /Cfci - fc^ functions. 
For the 2-point function this relation is the fluctuation-dissipation theorem of 
the last relation of Eq. (^7|) , while for higher n-point functions such a relation 
can be used to construct what have been advocated as generalizations of the 
fluctuation-dissipation theorem p. One can also see using Eq. (|^) that the 
n-point retarded functions in the two bases are related as 71raa---a = A^2ii---i- 
It is also of interest to examine other "mixed" functions, as for example 



'TIrraa-a{Pi,P2, ■■■,Pn) 

+ - coth''(a;p, / 2)/Ci2ii...i (Pi , P2 
1 



A^2211-l(Pl,P2, • • • ,Pn) + 
,Pn) + 



+ -COth''(Xp2/2)/C2111...l(Pl,P2, • • -jPn)- 



(70) 



Consider first the 3-point functions. Recent studies have examined the rela- 
tion between n-point functions of real-time formalisms and the analytically 
continued functions T^"'\zi, . . . , z„) of the imaginary-time formalism, where 
initially all energy components of Zi are in the discrete set of Matsubara 
frequencies on the imaginary axis. For the retarded products the relation is 

m 



|n|, 0, ig, |ii 



'TIrAA-a{,Pi,P2, ■■■,Pn) = ^211-l(Pl,P2, ■ 

= r(")(2;i +i{n-l)6,Z2-^ P2o 



■ ,Pn 

is, . 



P 



no 



ie).{n) 



Together with the complex conjugation relation of Eq. (46), this allows us to 
infer from Eq. ([70| ) the relation between the 3-point response function A^22i 
and various analytic continuations of the imaginary-time functions. Such 
a relation to the imaginary-time functions can also be worked out using 
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Eq. ( [70D for the 4-point response function /C22ii- For this we need, as well as 
Eq. ([7T|) , the relation between, for example, the mixed RA function TZaarr 
and various analytic continuations of the 4-point imaginary-time function 



la, 20 



We also note that Eq. ( ffOD can provide an interpretation of the "mixed" 
RA functions in terms of response functions in the Keldysh basis. As a result, 
certain causal properties (in time) can be inferred. For example, since the 
response relations of Eq. (|41|) are causal in nature, we conclude that each 
of the functions G'i222- -2(ii, ^2, ■ ■ ^2122-2(^1,^2, ■ ■ and G'ii22-2(ii, ^2, • • •) 
should only have ti or ^2 as the largest time 0, |T9[ . By a similar analysis for 



the non-amputated mixed functions Grraa - a as that leading to Eq. ([70|) , 
one could then deduce the corresponding causal properties for these mixed 
RA Green functions. 

We next consider the transformation between the FF functions and those 
of the RA basis. Writing this as 

(Pl, ...,Pn)= T^ihiPl) ■ ■ ■^a„/„(Pn)^/i.../„(Pl, • • • , Pn) (72) 

we find, for the choices a{p) = 1 and c{p) = e^n{x), d{p) = 77 in Eqs. 



51), the matrix 



Trf{p) Trf{p) \ ^ (I V . ^ ^ -1 
Taf{,p) Taf{.p) ) 




Oi-Po) [ I . (73) 



From this one sees that an F index corresponds to an R index for positive 
energies and to an A index for negative energies, and vice-versa for an F 
index. As well, with this relation one can translate the identities TZrr...r = 
= TZaa---a of the RA basis into the following set of relations in the FF 
basis: 

9{kQ) ■ ■ ■ 9{mo)9{-no) ■ ■ ■ 9{-ro)J^F -F P-pik, . . . ,m,n, . . . ,r) = 0, 
9{-ko) ■ ■ ■ 9{-mo)9{no) ■ ■ ■ 9{r^)J^F...F p...p{k, . . . , m, n, . . . , r) = 0.(74) 

These relations have a direct causality interpretion: if an F index represents 
a particle (po > 0) travelling forward in time and an antiparticle (po < 0) 
travelling backwards in time, and an F index represents the time-reversed 
situation, then Eq. ([TiD, together with the complex conjugation relation of 
Eq. (|55|), states that one cannot have particles and antiparticles either all 
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disappearing into the vacuum nor all appearing out of the vacuum. Such an 
interpretation also holds for arbitrary choices of the parameters a{p), c{p), 
and d{p). 

With the inverse of the transformation of Eq. ( ffSD we can, given relations 
such as Eq. ([7T|) relating the RA functions to the analytically continued 
imaginary-time functions T^"'\ examine how the FF functions are related 
to the imaginary-time functions. To simplify the notation, from now on we 
omit from T^"'\zi, . . . , Zn) the last complex argument Zn, which by energy- 
momentum conservation is minus the sum of the others. We first consider, 
with p + q = 0, the 2-point function J^ff{p, q) = —^pp{p, <?), for which we 
find 

:Fff{p, q) = ^(Po)r(^)(p; q) + 9i-po)T^^\p; g) = 

= e{po)r^^\z ^ po + ^e) + e{-po)r^^\z ^ po - ^e) = r^^\z po + ^epo), 

(75) 

as expected 0. The other two 2-point functions J^ff{pa) ~ ^pp^P^l) 
vanish identically, as can be shown by multiplying Tppip, g) by 1 = 6'(po) + 
6{—Pq) and then using the causality relations of Eq. (|73p. 

We next consider the 3-point function. With p + q + r = 0, we first 
examine J^fff{p, q, r), which using Eqs. PB| , [7I| , |T3D can be written as 



J^FFF{p,q,r) = 9{-po)9{-qo)T'^^\zi ^ po - ie,Z2 ^ go - i^) 
+0i-po)0{-ro)r^^\zi ^po- le, Z2 ^ qo + 2ie) 
+e{-qo)e{-ro)T^^\zi po + 2ie, Z2 ^ qo - le) 

n[—Xr) 

+^(po)^(ro)%^^r(3)(;^i po + ^e, ^2 - go - 2^e) 

+^(go)^(ro)^%^r(=^)(^, ^po- 2.e, - go + ^e). (76) 
n{-Xp) 

In this form Tfff appears to be a linear combination of the three possible 
analytic continuations of the 3-point imaginary-time function and their com- 
plex conjugates, but the presence of the 9 functions allows the replacement 
in Eq. ( |76D of all the V^^\zi,Z2) functions by a single analytically continued 
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function, T^^\zi ^ po + iepo, Z2 ^ qo + isqo) 0, 0]- Up to a normalization 
factor, then, this result could have been anticipated from the 2-point func- 
tion J^pp{p, q) of Eq. (|75D. In fact, a similar analysis for the n-point function 
J-'f-f{Pi, ■ ■ ■ ,Pn) indicates that this "Feynman" function in general can be 
written in terms of a single analytically continued imaginary-time function 
with ^ + iepiQ. 

One can also analyze the "mixed" FF functions in a similar way. For 
example, for the 3-point function ^PFpiPAi''^) "we find 

^FFpiP, Q, r) = 9{qo)9{ro)T^'''^ {zi ^ Po - ie, Z2 qo + 2ie) 
+e{po)e{ro)T'^^\zi ^po + 2ie, Z2 ^ qo - te) 

n[-Xq) 

-^^(-Po)^^(-ro)^^%^r(=^)(zi ^ po - 2re, z^ ^ q, + is). (77) 
n[-Xp) 

However, unlike the "pure" Feynman n-point function J^ff---Fi there does 
not appear to be any simple interpretation of these mixed FF functions in 
terms of single analytic continuations of the imaginary-time functions. 

Let us compare the number of independent functions present in all three 
bases. Taking into account the complex conjugation relation of Eq. (^) 
in the RA basis and Eq. (|55|) in the FF basis, as well as Eq. ( pT]) for the 
Keldysh basis, all three bases are seen to contain one independent 2-point 
function, which can be taken as the retarded function TZra = A^2i or the 
Feynman function TpF- For the 3-point function the RA basis has three 
independent functions, which could be taken as the three retarded functions 
T^RAA, T^ARA, and TZaar- The Keldysh basis contains these three retarded 
functions - /C211, /C121, and /C112 - but also four additional functions: /C222, 
A^i225 ^^^2125 and /C221- However, relations such as Eq. ( [fO|) and that following 
from TZrrr = in Eq. ( pD| ) shows that these latter four can be written as 
linear combinations of the retarded functions and their complex conjugates. 
Similarly, taking into account complex conjugation, the FF basis appears to 
have four 3-point functions, which could be taken as J^fff, ^fffi ^fff-i 
and J-'fff- However, relations such as Eqs. (^, |7^ show that all four can be 



written in terms of the three retarded products and their complex conjugates. 
Alternatively, recalling that these "mixed" FF functions are constrained by 
the relations of Eq. ( [7^ ) forcing energy to flow in only certain directions. 
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one could in this context count as an independent function each distinct 
orientation of the flow of external energy. In this way the RA basis has 
3 X 6 = 18 such independent oriented functions, while the FF basis also has 
1x6 + 3x4 = 18 independent oriented functions. One thus concludes that 
all three bases contain the same number of independent 3-point functions. 

Various examples could serve to compare the convenience of these bases 
in practical calculations. A simple and illustrative such example is the one- 
loop self-energy diagram of Eq. p7| ) with a bare +/— coupling constant 
g of Eq. (^5|) independent of momentum. The corresponding example in 
the imaginary-time formalism using the retarded analytic continuation z — >■ 
Po + ie has been explored by Weldon |23|. Let us then find from Eq. (^) 



the retarded self-energy. We first consider the Keldysh and RA bases, for 
which we calculate /C2i(p, —p) = T^RAip, —p)- In the Keldysh basis the non- 
vanishing bare vertices needed are ^211 and permutations of its indices, since 
D22{p) = 0, and so of all the possible terms in Eq. (0) only two survive. 
In the RA basis the relevant non-vanishing bare vertices are Qraa and qarr 
and permutations of their indices, since Drr{p) = = Daa{p)- Three 
possible contributions then remain in Eq. (|27|), but two of these disappear 
upon integration over because of the presence of terms like /S.R{k) /S.R[q) 
or A^(fc)A^(g) with poles in the same half-plane |jl5|. In both bases one 
thus obtains, with q = k — p and neglecting any 7-matrix structure. 



^21 (p, ~p) = TlRAip, -p) = 



(2vr) 



+ ^A{q) 



[Aniq)-AAiq)] + 
[An{k) - AA{k)]] , (78) 
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which agrees with previous results [|1^, | 

In this calculation the choice of free parameters b{p) = 1 or l/\/2 in 
the Keldysh basis and a{p) = 1 or —r]n{x) in the RA basis are easiest to 
use. The form of the result of Eq. ( |78[ ) is convenient for doing the k^ in- 
tegration, since for free propagators Aji{k) — AA^k) = 2ne{kQ)6{k'^ — w?). 
Such a form also arises for higher n-point functions at the one-loop level 
and, with suitable modifications, also at higher loop orders [15]. One might 
have anticipated that the RA basis would be significantly easier to use com- 
pared to the Keldysh basis: there are only two RA propagators Ar and 
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A^, whereas the Keldysh basis also contains the dependent combination 
As ~ coth''(x/2)[Ai^ — A^], and also two RA n-point functions TZrr.-.r and 
T^AA-'A vanish identically in general, while in the Keldysh basis only /C11...1 
does so. These features of the RA basis are certainly an advantage at higher 
loop orders. However, the vanishing of the bare Keldysh vertices with an 
even number of "2" indices, together with the particular form of A5 which 
arises directly in the result of Eq. (^), simplifies calculations in the Keldysh 
basis considerably. At least to lower orders, then, the use of the Keldysh 
basis and the RA basis as outlined here are comparable in difficulty. 

In the FF basis the retarded self-energy could be found as 9{po)J^ff{p, —p)- 
To show the final result of this calculation for Eq. (|27|) has the form as 
Eq. ([78D is somewhat tedious. An easier problem is to compute the imagi- 



nary part of this self-energy, which is discussed in detail by Weldon in the 
imaginary-time formalism . When this is found using Eq. ( P7D in the FF 
basis, terms such as 

J dko f{ko,qo) [AFik)AF{q) + AUk)AUq)] = 
= J dko /(go,ro) [A+(A;)A-(g) + A-(A;)A+(g)] , 

I dko g{ko, qo) [AFik)A*M + A^(fc) A^(g)] = 

= J dko 9iko,qo) [A+(A;)A+(g) + A-(fc)A-(g)] , (79) 

arise [|l7|, where, for example, A^{p) = 2'k9{±po)S{p'^ — w?) for bosons. 
Thus, 9 functions of energy automatically appear which, when applied to 
the vertices gfgh{p, q, r), lead to simplifications due to the normalization and 
complex conjugation relations described previously for the FF functions. 
The result in the end, of course, agrees with the corresponding one of Weldon 



and others |T|, |2|, |4| 



For this particular calculation of the imaginary part the choice of param- 
eters c{p) = e^n{x) and d{p) = t] in the FF basis is somewhat easier to 
use. However, in this example other methods such as the imaginary-time 
formalism or the RA or Keldysh bases are simpler, although this may not 
be so in more complex situations. In this regard, we remark that the steps 
used in this calculation are suggestive of the Cutkosky or cutting rules for the 
imaginary parts of graphs as derived in the approach of 't Hooft and Veltman 



[O. These rules at finite temperature using the +/— Dab propagators have 
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been discussed in Refs. [13, |25|, p6l ETl, and it is straightforward to derive 



the corresponding rules for the FF basis. A feature worth noting about such 
a derivation is that, because the zero temperature Feynman propagator and 
its complex conjugate are involved, the "non-cuttable" graphs that arise be- 
yond the one-loop level at finite temperature when the Dab propagators are 
used will not be present in the FF rules; for this the "causality" relations of 
Eq. (|7^ play a crucial role. 



7 Discussion 

In many instances physical quantities are more naturally expressed in terms 
of Green functions other than the usual time-ordered products. At finite 
temperature it is important to be aware of which function is needed, as dif- 
ferences between the various types of products is more pronounced than at 



zero temperature [0, Although in principle it is possible to extract all 
such information from the original Dab basis of the real-time formalism, in 
practice this may not be the most convenient method, especially in light of 
the constraints of Eqs. ^2|) among the components of Dab- For this reason 
formalisms such as the Keldysh and Retarded/Advanced bases have been con- 
structed which involve explicitly, in particular, the retarded propagator. The 
transformation of Eq. (Q) from the Dab basis to a generic transformed Dxy 
basis provides a unifying framework for these various real-time formalisms, 
in that the difference between them can be attributed to the different ways 



the constraints of Eqs. (|TT], |T2D are implemented. 



We have studied here three classes of such transformations - the Keldysh 
basis, the retarded/advanced RA basis, and a Feynman-like FF basis - which 



in a sense are attempts to implement the constraints of Eqs. ([11], "eco- 
nomically" . By means of the transformation between the bases we have also 
seen how information in one basis can be translated into the equivalent in- 
formation in another basis; this enabled us, for example, to examine the 
relation between various real-time functions and certain analytically contin- 
ued imaginary-time functions. Finally, we compared some aspects of these 
bases as they arise in practice. Although some of the relations in the FF 
basis such as complex conjugation can be chosen to be simpler than the cor- 
responding ones in the RA basis, it appears in general that the conditions 
of Eqs. ([Tl|, 0) are implemented more efficiently in actual calculations in 
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the RA basis ITR |28| 



30 



or even, comparable at least to lower orders, 
in the Keldysh basis. Thus, for practical purposes, it appears that of these 
three bases the RA or perhaps the Keldysh basis is easier to use in general, 
although this could depend on the particular apphcation and may also be a 
matter of taste. 
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